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Space-Constrained Massive MIMO: Hitting the 
Wall of Favorable Propagation 
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Abstract —The recent development of the massive multiple- 
input multiple-output (MIMO) paradigm, has been extensively 
based on the pursuit of favorable propagation : in the asymptotic 
limit, the channel vectors become nearly orthogonal and inter¬ 
user interference tends to zero CD- In this context, previous stud¬ 
ies have considered fixed inter-antenna distance, which implies an 
increasing array aperture as the number of elements increases. 
Here, we focus on a practical, space-constrained topology, where 
an increase in the number of antenna elements in a fixed 
total space imposes an inversely proportional decrease in the 
inter-antenna distance. Our analysis shows that, contrary to 
existing studies, inter-user interference does not vanish in the 
massive MIMO regime, thereby creating a saturation effect on 
the achievable rate. 

Index Terms —Antenna arrays, favorable propagation, massive 
MIMO. 

I. Introduction 

Massive multiple-input-multiple-output (MIMO) systems 
with hundreds of antennas are envisaged for the next genera¬ 
tion of base stations (BSs) to meet the ever increasing quality 
of service (QoS) demands in a power efficient manner [Q- 
E). Existing studies have shown that a key characteristic of 
massive MIMO is that the inner product of two distinct channel 
vectors between the BS and the mobile users tends to zero as 
the number of BS antennas increases ID, 0. In other words, 
the multi-user spatial streams become asymptotically pairwise 
orthogonal. Under this condition, coined favorable propa¬ 
gation, low-complexity precoding/detection (e.g. maximum- 
ratio-transmission (MRT)/maximum-ratio-combining (MRC)), 
can achieve near-optimal performance (T|. 

A major obstacle towards increasing the numbers of anten¬ 
nas in practical systems is the limited available physical space 
in both BSs and mobile devices. The placement of antenna 
elements at more than half a wavelength apart is, in general, 
considered to secure minimal correlation between the commu¬ 
nication channels, allowing the exploitation of the full system 
spatial diversity. On the other hand, a dense deployment of 
antenna elements increases the spatial correlation, thereby 
deteriorating the system performance (4|, (5|. Accordingly, 
for a fixed total physical space, this introduces a tradeoff in 
the spatial diversity. With an increasing number of antenna 
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elements, along with the increasing signal sources, spatial 
correlation increases as well, thereby inducing a saturation 
in the spatial diversity benefits. This observation motivates 
a study on whether the favorable propagation characteristics 
hold for the case of space-constrained massive arrays. 

In line with the above, in this paper we study the sta¬ 
tistical properties of the inner product of channel vectors 
for a system equipped with a massive uniform linear array 
(ULA) of fixed total space, where the inter-antenna distance 
is inversely proportional to the number of elements and the 
antenna correlation is directly dependent on the number of 
antennas. In contrast to a, a, where the sum rates were 
derived for particular precoding schemes, here we elaborate 
on the fundamental statistical properties of the massive MIMO 
channel for a space-constrained antenna deployment. Our 
analysis shows that the inner product of two distinct channel 
vectors, which intimately represents inter-user interference, 
converges to a non-zero value; this causes a saturation of the 
achievable sum rates under MRT with an increasing numbers 
of antennas. In fact, it is analytically shown that the first and 
second-order moments of this inner product converge to non¬ 
vanishing values which depend on the total physical space 
of the BS array. Note that the findings under this scenario 
stand in fundamental contrast with those of (T], (3), (51. which 
demonstrated that inter-user interference converges smoothly 
to zero for physically unconstrained antenna arrays. 


II. System Model and Achievable Sum Rate 
A. Channel Model 

Consider a downlink multiuser MIMO system with a N- 
antenna BS and K single antenna receivers with N > K. 
The physical space for the antenna array at the transmitter is 
constrained to do A, where A is the carrier wavelength, and do 
is the arbitrary physical constraint in units of wavelengths. The 
channel is assumed to experience frequency flat fading, such 
that the received signal at the k -th terminal reads as 

K 

y k = x/phf WfcSfc + y/p ^2 h fw jSj + n k (1) 

where (-) H denotes Hermitian transpose, n k ~ CA/”(0,1) is the 
additive white Gaussian noise vector, p is the average signal 
to noise ratio (SNR), represents the downlink channel to 
the fc-th terminal, w/. and s k are the unit-norm beamforming 
vector and the information symbol from the BS for the k -th 
terminal, respectively. For a fixed total space of do at the BS, 
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the inter-antenna separation is given a 


H 


d = 



( 2 ) 


The generic channel model using steering vector representation 
can be expressed as 


L 

h k = /3k,og{0 kfi ) + ^2 Pk,ig{0k,e) (3) 

e=i 


where the first term models the line-of-sight (LOS) component 
with path gain (3k,o and the second term accounts for the L 
multipath components with path gains f3k,e.,(- £ [1,L]. The 
steering vectors gk are modeled as 


g (Ok,e) = 


l,e 


-i2n4- sin 6, 


k,i 


0 -i2n(N-l)% sin 8 ktl 


1 T 


(4) 


where (-) T denotes the matrix transpose, i = y/—l and 0 k 
is the angle of departure (AoD) for the fc-th terminal. In line 
with El. El we focus on the uniform random LOS channel 
with /3k, o = 1 and f3 k ,i = 0 for all t ^ 0, for which we have 


h fc = 


l,e“ 


V2-k{N sin 9 k 


(5) 


It can be seen that now the channel response and the resulting 
antenna correlation are a function of the number of antennas 
N. We assume that all terminals are randomly distributed 
within a circle-shaped cell with radius R. Following 0, 0, 
sin 9k, k = 1 are assumed to be uniformly distributed 

within the interval [—1,1]. 


B. Achievable Sum Rate 

To keep our analysis insightful, we assume that both the BS 
and users have perfect channel state information. According 
to 0, the achievable ergodic rate of the /c-th user can be 
expressed as 


Rk =£{log 2 (l + SINR fc )} (6) 

where E{ ■} denotes the expectation operator and 


SINR fc = 


Pl h f w fcl 2 

l + PEf= lj¥ Jhfw,f 


(7) 


Considering MRT where the precoding vector is given as 
Wfc = h k/y/N, the achievable ergodic rate of the /c-th user 
can be rewritten as 


R k =E{ log 2 1 + 


> log 2 1 + 


pN 


l + ^Ef=i, j¥fe |hfh,l^ 
1 


( 8 ) 




= Rl 


(9) 


^ote that, depending on the specific ULA configuration, the total space 
may be split in N — 1 equal separations, giving d = A. For notational 
convenience we herein consider N such separations, assuming d/2 additional 
space at each end of the array, which for large N will give a negligible 
difference between the two models. 


where © has been obtained by applying Jensen’s inequality 
on the convex function of the form log 2 (l + 2). The above 
expression illustrates that the inner product, h].' h, ; , between 
two distinct channel vectors determines the resulting inter-user 
interference and, consequently, the achievable rate of each 
terminal. In the following, we will investigate the statistical 
properties of this inner product and show that, in the asymp¬ 
totic limit as N —> oo, its first two moments converge to 
non-zero values. 


III. Non-Favorable Propagation Conditions 

In this section, we derive exact and asymptotic expressions 
for the first and second-order moments of the inner product 
h^hj, which determine the inter-user interference power. We 
begin with the exact analysis, which holds for any finite N, 
and present the following proposition: 

Proposition 1. For physically constrained ULAs with a total 
space of do in the LOS channel, the first two moments of the 
inner product hj^hj are equal to 


N-l 

.Ejh^lx, j = ^2 sine 2 (am) 

771=0 



N-l N-l 

= ^ 2 ^2 s ^ nc2 ( a ( m i ~ m 2 )) 

?71i=0 7712 —0 


/N-l 


- E 


sme (am 


V, 771=0 


( 10 ) 


(ID 


where a = 27t^ for brevity and var(-) returns the variance 
of a random variable. 

Proof: Let us define u k = sin 6^. Then, we have 


N-l 

#{hf hj} ='% 2 e J (12) 

771=0 


from which, for a uniformly distributed u k in the interval 
[—1,1], it can be seen that 


sin(27r^m) 

2t r^m 

The desired results in (ITOt -dUl) follow trivially after some 
basic algebraic manipulations and noting that 



var 


jhfhjj =£ j|hfh,-| 2 } — E{ hfh,-} 


(13) 


We now turn our attention to the massive MIMO regime, 
by letting N —> 00 . Under these conditions, we obtain some 
very important insights. 

Proposition 2. For physically constrained ULAs with a total 
space of d 0 in the LOS channel, the first two scaled moments 
of the inner product li / // li ; and the square scaled inner product 
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Fig. 1. Mean channel vector inner product E { -^hj^h j } vs. N, for massive 
MIMO with limited/unlimited physical space. 



Fig. 2. a) Mean channel vector inner product E { jrh^hj } vs. do and b) 
variance of channel vector inner product var hj} vs. do, for massive 

MIMO with limited/unlimited physical space (N = 200). 


|hfh,f tend asymptotically (TV —> oo) to 

¥5 |h " h ' |2 }^i- £ 

var/ 1 —y —— fl— 77— — 77 —^ 

\ N j 2 cIq \ 2N 8do J 

1 1 v^°° sin 2 (am) • . • , 

where e = ?,•> > ™ 1 --—- is a correction term. 

4-7HaQ m 

Proof: See Appendix HT1 ■ 

The above results indicate that the channel vectors do not 
become pairwise orthogonal in the massive MIMO limit. 
Instead, their mean inner product converges to a non-zero 
value that is inversely proportional to the total physical space 
of the array. In other words, the smaller the physical length of 
the BS array, the more we deviate from favorable propagation 
conditions. Yet, if do grows analogously with TV, the mean 
inner product will converge smoothly to zero, as predicted in 

m, 0,0. 

From IS Eq. (18)—(19)], we know that for the case of 
unlimited physical array space, 

(17) 

var {^ h " h J } = ^-jT (18) 

Comparing (IT4l) and (fl 6 ]> with ( 1 1 7| i-( l I X| i. we can infer that the 
corresponding expressions coincide if do = 4 - in which case 
both moments will converge smoothly to zero. On the contrary, 
if the physical array space is fixed (i.e. it does not grow with 
TV), the mean and variance of the inner product, which indicate 
the level of inter-user interference, will converge to non-zero 
constant values that are inversely proportional to do- 



(14) 

(15) 

1 

4TV^~ 6 

(16) 


IV. Numerical Results 

This section illustrates the analytical performance and re¬ 
sults of Monte-Carlo simulations for the massive MIMO 
scenario under study in comparison with the case of massive 
MIMO with unlimited physical space. The simulations assume 
a linear antenna array at the BS, I\ = 10 users and the steering 
vector channel model of ([5]). 

Figure [I] shows the simulated and analytical mean inner 
product (flOt of the channel vectors for an increasing number 
of antennas TV, along with the theoretical limit of (fl4l >. The 
cases for do = 4 and do = 10 are shown and it can be seen that 
m yields an exact match to simulation. More importantly, 
it is evident that, unlike the case of unlimited physical space, 
for the space-constrained deployment the inner product does 
not converge to zero, but to -j4- as shown in our analysis. 
This implies that interference is not completely nulled, which 
makes MRT precoding non-optimal for this scenario. 

In Fig. [3a), the simulated and analytical mean inner product 
of the channel vectors is shown for an increasing total ULA 
space do- The mean inner product is also shown for the ULA 
deployment without space constraints for reference. A non- 
negligible deviation between the two models can be seen, 
which becomes significant as the total space d 0 decreases. 

In Fig. |2|b), the simulated and asymptotic variance of the 
inner product of the channel vectors is shown for am increasing 
total ULA space do. The analytical asymptotic limit illustrates 
expression (fl6l) . which yields excellent tightness with the 
simulation results even for small do- Again, the case of ULA 
deployment without space constraints is shown for reference 
and a deviation between the two models can be seen. 

Finally, Fig. [3] illustrates the effect of the above observations 
on the achievable sum rate of the system under MRT, for a 
transmit SNR of p = lOdB where we have accounted for 
additional shadowing and propagation losses of 30dB. The 
Jensen’s lower bound in <(9]) can be easily evaluated in closed- 
form using the results of Proposition 1; unfortunately, /T/ is 
relatively loose for the channel model considered herein, since 
var j | 2 } is much higher compared, for example, to 

the Rayleigh fading case considered in 0 . The achievable 
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Fig. 3. Achievable sum rate for massive MIMO with limited/unlimited 
physical space. 


sum rates are shown for the cases do = 4,10 and compared to 
the deployment over unlimited physical space. Contrary to the 
latter case, a saturation in the sum rates can be observed for the 
space-constrained deployment, where, as expected, the ceiling 
of the achievable sum rate starts diminishing with increasing 
total space do- 


Conclusion 

The favorable propagation characteristics of massive MIMO 
have been examined under space-constrained antenna deploy¬ 
ments. Contrary to current understanding, it was shown that 
when increasing the numbers of antennas that are deployed in a 
fixed total space, the inner product of the channel vectors tends 
to a non-zero value. This phenomenon implies that inter-user 
interference does not vanish in the massive MIMO regime, 
thereby making MRT non-optimal. The above discussion poses 
a fundamental divergence from existing intuition and motivates 
a further study of the resulting space-constrained deployments 
and the development of suitable transmission schemes. 


Appendix II 

Proof of Proposition [2] 

The proof of ( 11 4b follows trivially by recalling (1 1 01 1 and 
Lemma 1, where the condition on a holds for large N. 
Regarding the square inner product E ||h^hj|“|, this can 
be expanded as follows: 


N-l N-l 


E 


{| h f h jf}= E E sine 2 (a(mi - m 2 )) 


772 1=0 7722=0 

N-l N-l 


= »+Y Y 


m 1—0 r?22 = 0 

7722 / 7721 


1 — cos( 2 a(mi — m 2 )) 
2a 2 (mi — m^) 2 


N-l N-l 


= N + 


sE E 


2 a 2 ^ ^ (mi - m 2 ) 2 

772l=0 772 2 = 0 

7722 # 7721 


N-l N-l 

Y Y 


cos( 2 a(mi — m 2 )) 


1 

2 a 2 2 -—' 2 -—' (mi — m 2) 2 

m 1=0 m 2 = 0 

7722 # 7721 


( 21 ) 


After some tedious but straightforward manipulations, it can 
be shown that, for N —> 00 , the first sum in (fTi~l) converges to 


N-l N-l 

Y Y 

m i=0 m 2 = 0 
m 2 5 ^ mi 


(mi — m 2) 2 6 


2 00 1 

,7T \ - 1 




k=l 


Likewise, the second sum in (fTil converges to 


N-l 


N-l 


Y Y 

mi=0 m 2 = 0 

m 2 7 ^ mi 


cos( 2 a(mi — m 2 )) 
(m 1 — m 2) 2 


^2iVp-,ra + L)-£t!2<M. (23) 


fc =1 


Combining (f2l1) — (f23l) and simplifying, yields the desired re¬ 
sult. Finally, the asymptotic result for the variance of the inner 
product is a direct combination of (ITll) . (IT5l) with (IT3l) . 


Appendix I 
Useful Lemma 

Lemma 1. Defining sinc(x) = sin(x)/x and for 0 < /? < 7 r, 
we have that 

Y sinc2 = \ ( X + • (19) 

Proof: Using standard trigonometric properties, the infi¬ 
nite summation in ( IT9l ) can be expanded as follows: 


lxj - LXJ 

Y sinc2 = 1 + ^2 Y 

772=0 772=1 


1 


E 


cos(2/3m) 
2 /3 2 m 2 


( 20 ) 


The first infinite sum in ( l20t converges to 7 t 2 /6 while the 
second one can be evaluated with the aid of GjEq. (5.4.2.12)]. 
Putting everything together, yields the desired result. ■ 
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